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Abstract 

There are several physically motivated density matrix functionals in the literature, 
built from the knowledge of the natural orbitals and the occupation numbers of the one- 
body reduced density matrix. With the help of the equivalent phase-space formalism, 
we thoroughly test some of the most popular of those functionals on a completely 
solvable model. 

1 Introduction 

At the Colorado conference on Molecular Quantum Mechanics (1959), Coulson pointed out 
that in the standard approximation the two-body density matrix 72 carries all necessary 
information required for calculating the quantum properties of atoms and molecules [1]. 
Because electrons interact pairwise, the main idea consists in systematically replacing the 
quantum wave function by the two-body reduced matrix (a function of four spatial variables), 
which may be obtained by integration of the original A-body density matrix (a function of 
2 A spatial variables). The A-representability problem for the two-body reduced density 
matrix has proved to be a major challenge for theoretical quantum chemistry [2]. 

In the last fifteen years there has been a considerable amount of work on Ansatze for 
the two-body matrix in terms of the one-body density matrix 71. Starting with the pioneer 
work by Miiller [3], rediscovered in [4], several competing functionals have been designed, 



partly out of theoretical prejudice, partly with the aim of improving predictions for particular 
systems: among others, the total energy of molecular dissociation [5] , the correlation energy 
of the homogeneous electron gas [6] and the band gap behavior of some semiconductors [7] . 

Two-electron systems are special in the sense that 72 can be reconstructed "almost ex- 
actly" in terms of 71. Namely, let us express 71 by means of the spectral theorem in terms 
of its natural spin orbitals {4>i{r)} and its occupation numbers {n^}. The ground state of 
this system (which is of closed-shell type) admits a one-density matrix: 
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x,x') = (titi' + Wi')7i(»-,r') = ( Tl 1V X>&(r)#(r'), (K 



where x = (r, s) stands for the spatial and spin coordinates. The natural occupation numbers 
satisfy Yli n i = 1> w ^ n < n « < 1- Mathematically this is a mixed state. The original two- 
density matrix is given by the Shull-Lowdin-Kutzelnigg (SLK) formula [8,9]: 

l2 (x 1 ,x 2 ;x' 1 ,x' 2 ) = (til 2 - I1T2) (ti4 2 ' ~ IM^^ Mri)Mr2)<f>*(r[)<f)*(r 2 ), 

with coefficients q = ±y/n~i . (2) 



Though the expression is exact, the signs of the q still need to be determined to find the 
ground state. Note that 

71(33; x) = 2 J -f 2 (x,x 2 ;x',x 2 ) dx 2 . (3) 

This condition becomes a sum rule which, as we shall see, may be satisfied or not by proposed 
density matrix functionals. Among other conditions, 72 is Hermitian: 72(^1, x 2 ; x[, x' 2 ) = 
j2(x[,x 2 ;Xi,x 2 ), and antisymmetric in each pair of subindices: 

7 2 (cci, x 2 ; x[, x 2 ) = -72(^2, Xi; x[, x' 2 ) = -72(^1, x 2\ x 2 , x [)- 

In the early years of the theory, Heisenberg invented an exactly solvable model, here 
called harmonium, as a proxy for the spectral problem of two-electron atoms [10]. It exhibits 
two fermions interacting with an external harmonic potential and repelling each other by a 
Hooke-type force; its Hamiltonian, in Hartree-like units, is 

rj \Pl\ 2 , Irf . k 2 2 5 2 

H = + + ^(\r!\ + \r 2 \ ) - -r 12 , (4) 

where 5 > and r 12 := |ti — r 2 \. Many years later, Moshinsky [11] came back to it with the 
purpose of calibrating correlation energy — see also [12-14]. Also, Srednicki [15] used the 
harmonium model to study the black hole entropy, proving its proportionality to the black 
hole area. 

Recently, within the context of a phase-space density functional theory [16], here called 
WDFT, the alternating choice of signs in (2) has been shown to be the correct one for 
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the harmonium ground state [17-19]. Mathematically, density matrix functional theory 
(DMFT) and density functional theory on phase space are equivalent: see the next section. 
Thus, its density matrix functional (2) is nowadays known exactly. Some of lower excited 
configurations of harmonium are also of much current interest [20-22] . 

The integrability and solvability of harmonium enables one to test accurately how pro- 
posed density matrix functional behave for this particular system [23]. There is another 
particularity of harmonium: for its ground state, the Miiller functional, evaluated on the 
exact one-body reduced density matrix, yields the correct value of the energy [24]. Here we 
confirm by a different method this surprising coincidence, and we catalogue the predictions 
for harmonium by several proposed two-body functionals, measured against the exact model. 

In section 2 we briefly recall the analytical phase-space treatment for the harmonium 
ground state; this also helps to introduce the notation. In section 3 the Hartree-Fock 
and Miiller functionals of the one-body density are discussed in the context of harmonium. 
Section 4 starts the systematic comparison of several other approximate functionals proposed 
in the literature; we compute the error in the interelectronic energy value given by each of 
them for the exact family of ground states parametrized by (k, 8). They all behave worse 
than Muller's. We append some concluding remarks. 

2 Wigner natural orbitals for the harmonium ground state 

The basic object of WDFT for a two-electron atom is the one-body quasiprobability d±. Let 
us look at the two-body quasiprobability d 2 . Given any interference operator | X I / )( < I ) | acting 
on the Hilbert space of the two-electron system, we denote 

r 2 ; Pi,P2] ft, ft; ft', ft') (5) 



These are 4x4 matrices on spin space. When the interference operator corresponds to a pure 
state = <3>) we speak of Wigner quasiprobabilities. In this case, the functions (5) are real, 
and we write d 2 for P 2 . Its integral equals 1. The extension of this definition to mixed states 
is immediate. The corresponding reduced one-body functions are found by integration: 



On spin space these are 2x2 matrices. When \I> = <3> we write d\ for P\. The integral of 
this quantity equals 2. The associated spinless quantities are obtained by tracing on the spin 
variables: 




*(n - z u r 2 - z 2 - qi , ft ) + z l} r 2 + z 2 ; e *(pi-*i+w«) dZl dz2 . 



A*$(n;pi;ft;ft') = 2 / P 2 **(»*i, r 2 ;pi,p 2 ; ft,ft;ft', ft) dr 2 dp 2 dq 2 . 
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The marginals of d 2 give the pairs densities P2(Ti,r 2 ), ^(Pijf^)- The marginals of d\ give 
the electronic density, namely p(r) = J di(r,p) dp, and the momentum density 7r(p) = 
J d\(r,p) dr. It should be obvious how to extend the definitions to iV-electron systems and 
their reduced quantities; the combinatorial factor for d^ >->■ d n is (^). 

Putting together the equations (2) and (1) with (5), one arrives [17] at: 

d2{.r l ,r 2 ;puP2^i^2^v^2') = (spin factor) x ^ ^ Xij{ r u Pi)Xij ( r 2; P2), (6) 

and di(ri;pi; <?i, ft') = 2 ^ d 2 (T*i, r 2 ;pi,P2; ft, ft; ft', ft) dq 2 dr 2 dp 2 

= (titr + iiiv) Yl Ui Xi(ri;pi), 

i 

where n« are the occupation numbers with < rii < 1 and = 1> the X«j the natural 

Wigner interferences and x« := X«« denote the natural Wigner orbitals; the spin factor is 
that of (2). Evidently (titi' + iiiv) is a rotational scalar; thus we replace it by 2 in what 
follows. 

Introducing extracule and intracule coordinates, respectively given by 

R= -^=(7*1 + r 2 ), r = -j=(ri-r 2 ), 

P = -^(Pi+P2), P=^(Pi-P2), 

the harmonium Hamiltonian (4) is rewritten: 

rr rr rr P 2 ^ 2 B 2 p 2 A 2 

where u; := and [i := V&; — 5. We assume 5 < k. The energy spectrum for harmonium 
is obviously (N + §)w + (N + |)/x and the energy of the ground state is E = |(w + //). For 
this configuration, the (spinless) Wigner two-body quasiprobability is readily found [16]: 

At ^ 1 ( 2Hr \ ( 2H A M 

d2{r 1 ,r 2 ;p 1 ,p 2 ) = — exp exp . (7) 



The reduced one-body phase-space quasiprobability for the ground state is thus obtained: 



d 1 (r;p)=2 J d 2 (r,r 2 ;p,p 2 ) dr 2 dp 2 = ^ 



2 / 4w/i \ 3/2 _ 2r . 



3 v (u + Li) 2 



e -2r 2 a; M /(w+ M ) e -2p 2 /(^+M) i 



Its natural orbital expansion, with i integer > and the corresponding Laguerre polyno- 
mial, reads [17]: 

.2 _ „ . 4 V^ / V 7 ^ - VJ 1 \ 21 _. /1 ^ +2i 



c\ = ni = v - 9 -V ^= =:(l-t 2 )* > (8) 

Xi(n;Pi) =Xi(^i;pix)xi(yi;pij/)xi(^i;pu), where 

7T 
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Up to a phase, the functions x% determine the set of interferences: for j > k, 

Xik(x,p x ) = \ (-1)^ {2V^x 2 + 2pl/^) {j - k)/2 

x e'^-WLi- k (2^JJlx 2 + 2pl/ v fo^)e-^ x2 - p * / ^, (9) 

where d := &rctan(j) x /y/uJJIx). The L]^ k are associated Laguerre polynomials. The Xkj are 
complex conjugates of the Xjk- 

The SLK relation q = ±y/rTi must hold. There is the problem of determining the signs 
of this infinite set of square roots, so as to find the ground state. In principle, to recover d 2 
from d\ is no mean feat, since it involves going from a statistical mixture to a pure state. 
The advantage in the present case is that the result (7) is simple, particularly so on phase 
space, and known. With the alternating choice (unique up to a global sign): 

c l = (-r v /^ = v / r^Hy, (io) 

and the above fjk, formula (6) does reproduce (7). This was originally proved in [17] by 
organizing the series (6) in a square array and summing over subdiagonals; some special 
function identities come in handy at the end. 

Incidentally, new special function identities certainly lurk here: a natural idea in this 
context is to try to sum the SLK series differently. Consider for instance the sum on the 
first column: 

So := E^X^PrWr^) = ^ e'^^ £(-l)V±(4^)'/V fr <*+*> 

r r=0 



1 -t 



2 



e 



-(U}+Uj) 



7r 2 ' — ' r 

r=0 



oo 

Y^-^tu^e-^^y 



' * \-^ +u 'hxp{-2tU 1 U 2 e-^ + ^), 



n 2 



with the notations 

1/2 



Ui := [y/ujirf + p 2 / v /ujl] , ^ := arctan(p i / v /w/ir i ), 2 = 1,2. 

As recalled in [25], the exponential generating function of the Bell polynomials B n (y) is 
exp(y(e x - 1)) = Y^ =Q B n {y) x n /n\. The polynomials are given by B n (y) = YT m =o {1} V m > 
where {^} is the number of partitions of n into exactly m subsets. Therefore, 



exp(j/e x ) = e y ^B n {y) 

n=0 

and so 

f2 



X 



S Q = 1 -^ e -(^^) V B n (-2tU lU2 ) tW±±jW 
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We have not found a clear way ahead, however, for the summation of all columns. On the 
other hand, the true and tested minimization method works fine to derive (10), too [18]. 
Trivially, the same sign rule holds for natural orbitals of the garden variety (2). This is 
invoked in [26] without explanation. 



3 Hartree— Fock and Miiller functionals of the one-body density 

The total energy of an electronic system is of the form 

Etotal[dl, d 2 ] = -Skinful] + -E'extMl] + E ee [d 2 ], 

where the kinetic £kin an d potential E ext energies are known functionals of the spinless 
one-body Wigner reduced quasiprobability. In our case, recalling (8): 



Eki n [di] 
E ext [di] 



p 2 3w 3/i 3u 

— d l (r-p)drdp = — + — = — 



1 + 



1 - t\2 



ui 2 r 2 



, . . , , 3u 3u 2 
— di(r;p) drdp= — + — 



3co 

T 



1 + 



i + t 
1 + t\ 2 



The interelectronic repulsion energy E ee [d 2 ] is a functional of the two-body Wigner quasi- 
probability or indeed only of the pairs density: 



E ee [d 2 



d 2 (ri, r 2 ; P\,Pi)r\ 2 dr x dr 2 dp x dp 2 



5 f 

E ee [p 2 ] = -- / p2(r 1 ,r 2 )r 2 2 dr 1 dr 2 



3/i 3w 5 



3cu 8t(l + t 2 ) 



4 4/i 4 (1-t 2 ) 2 



We note that the kinetic energy of the system stays finite from /x = u (t = 0) to /i = 
(t — 1). The potential energies diverge as \x \. 0, in the strong repulsion regime; but their 
sum remains finite and equal to the kinetic energy, as prescribed by the virial theorem. 

In the language of this paper DMFT amounts to the search for functionals for p 2 in 
terms of d 1 , expressed through its Wigner natural orbitals and their occupation numbers. 
We have already indicated that the wide variety of functionals currently used in DMFT for 
computational purposes can be traced back to the functional proposed by Miiller [3]. Note 
first that, with an obvious notation, the exact phase-space functional for the present system 
is: 



(11) 



i,3>0 



^2n i p i (r l )p i (r 2 )+ ^ ^njTj XijMXijfa), 



i>0 



with pi(r) = Xi( r ) being the electronic density for the natural orbital Xi{ r iP)- This is 
correctly normalized by f p 2 (ri,r 2 ) dr±dr 2 = X^i n « = 1; m y i ew of f Xij( r )dr = when 
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% 7^ j. Translated into our language, the Miiller functional for the singlet is of the form 

pf(ri,r 2 ) =2j2n l Xt(r 1 )J2 n jXj(r2)-Y,\ m ^Xij(r 1 )xj l (r 2 ) (12) 

i j i,j 

= ^p(ri)p(r 2 ) - ^riiPiirJpifa) - ^ V n * n i Xij(ri) Xy(r 2 ). 

i ijtj 

We have used that Xij{ r ) = Xji{ r ) f° r re& l orbitals. More generally, Miiller considered nj 
with p + q = 1 instead of ^Jn \ n~j . Recently, the case \ < p = q < 1 has been studied [7]. 
The Miiller functional satisfies some nice properties; among them, the sum rule (3) and 
hermiticity. For Coulombian systems its energy functional is convex [27]. Nonetheless, 
antisymmetry fails. (We summarize properties fulfilled or infringed by each functional in 
Table 1.) 

Following Lieb [28], the Hartree-Fock approximation may be regarded as yet another 
functional of d 1 . This is given by 

pf F (ri,r 2 ) = -p(ri)p(r 2 ) - XX#( r 0#( r 2) -^mnjXijirJXijfa). (13) 

i i^j 

Expressions (11) and (13) coincide only when the occupation numbers are pinned to or 1. 
The cumulant p 2 — P^ F can be also computed [29]. The "best" Hartree-Fock state, in the 
sense of best approximation for the ground state energy with only one rii ^ 0, is given [19] 
by: 

7T b 

= -1 e -(R 2 +r 2 )y/(^+^)/2 e -(P 2 +P 2 )/y/(^+^)/2 
7T 6 

Use of the energy formulas for this state yields 

£ H f = 3vV + /i 2 )/2, 

and so the correlation energy is 

E c (ou,p) :=E -E HF = ^(w + p- v / 2(w 2 + /i 2 )) 
3 

so that E c (u,0) = -~{V2 - l)u. 

For small values of p, however, minimization by use of (13) gives lower values of the energy 
than i?HF [14]: the results by Lieb on the Hartree-Fock functional for arbitrary states of 
Coulombian systems do not apply here. However, it should be remembered that the sum 
rule fails for non-Hartree-Fock states. 
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4 Exact vs. approximate functionals for harmonium ground state 

Our strategy henceforth is simply to gauge the worth of the functionals by computing their 
respective values on the true ground state. As mentioned above, it has recently been found 
[24] that the Miiller functional yields precisely the correct energy values for harmonium 
when evaluated on the exact ground state — thus, for TV = 2, it is also overbinding for the 
harmonic repulsion just as for the Coulombian one [17,27], since the minimizing state for 
that functional will yield a lower value of the energy. 

Thus a feasible procedure is to compute the difference between the values given by the 
Miiller functional and each of the several functionals whose accuracy we want to study. We 
need only worry about the interelectronic repulsion energy; since all the relevant quantities 
factorize, for notational simplicity we shall work in dimension one. 



4.1 Miiller interelectronic energy 

From (12) and (9), we get: 



\fn r n s Xrs (n , Pi ) Xsr (r 2 , Vi ) 

r,s=0 

1 _m2,rr2\ t~ t .n/ —ill) s\ 



_i_ e -(uHu!) £ {2UiU2t y e -ue £ 4(2[/2) L , W) ^ 



l=r—s s>0 

oo 



7T 2 \ 1 — t 2 

l=-oo V 

= JL e -(^?+C/|)(l+t 2 )/(l-t 2 ) e 4UiU a t<x*0/Q.-1?) 
7T 2 

where = $i — $2; some well-known properties of Laguerre polynomials and modified Bessel 
functions have been invoked. In all, the spinless phase-space Miiller functional for the har- 
monium ground state is, using the notation U{ := (ri,pi), 

df( Ul ,u 2 ) (14) 

- 2 ( l ~ t2 \ p -(U?+Ul)(l-t 2 )/(l+t 2 ) _ ^ p -(U'i+U'i)(l+t 2 )/(l-t 2 ) 4U!U 2 t cos e/(i-t 2 ) 
" TT 2 \l+t 2 J TT 2 

In order to compute the interelectronic energy, we proceed with the mean value of the 
electronic separation: f r\ 2 d^(ui, u 2 ) du\ du 2 - For the first term in (14), we get: 

7T U + jl J UJ/JL 
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Figure 1: Exact vs. Miiller pairs density for harmonium at £ = 0, t = 0.3 and t = 0.5. The 
graphics show the dimensionless functions 7rp 2 /w (on the left) and np^f/u (on the right) in 
terms of ui 1 ^ 2 1\ and ui 1 ^ 2 r 2 . 



For the second term, we obtain: 



7T 2 



e -(Uf+Ui)(l+t 2 )/(l-t 2 ) e 4C/il/ 2 tcose/(l-t 2 ) d d 

cupl 



r 2 e-^i+r^+n) e nra(«-Ai) d r 2 dn = -- . 



12 1 

71 / UJ 



In the process we have obtained a sort of (spinless) "Miiller pairs density" for the true ground 
state, 



J ' " 7T UJ + fM 7T 
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Figure 2: Diagonal part 7rp^ 1 (r, r)/u and antidiagonal part irp^(r, —r)/u of the Miiller func- 
tional as functions of t and uj 1 / 2 r. 

whose predicted mean square value for the distance between the two electrons is 



Thus, both energies coincide: E ee = = — <5/4/i = (/i 2 — w 2 )/4yU. Note that the result is 
valid for any value of t. This is surprising because the shapes of p 2 and p^ 1 grow very distinct 
as t increases — see Figure 1. 

In summary, by a somewhat different method, we have confirmed the result of [24]. The 
coincidence does not hold for other values p,q ^ | in the Miiller approach. It may be 
considered fortuitous, because (15) and the exact pairs density (16) are rather dissimilar: 
for t > 0, the spinless two-body Miiller functional does not have a maximum at the origin 
in phase space, whereas the exact functional does. More precisely, as figures 1 and 2 show, 
the Miiller functional exhibits two maxima located at the antidiagonal sector of the density. 
Also, it sports negative values at some points. As pointed out in the original paper [3], this 
phenomenon is a consequence of the inequality ^/n] > rij satisfied by the natural occupation 
numbers of the system. Figures 1 and 2 display the negativity around the diagonal elements 
of the density. This indicates that the Miiller functional is also unphysical, in a subtler way 
than the Hartree-Fock functional [30]. 

4.2 Hartree— Fock interelectronic energy 

We use the following terms, computed in [19]: 




The same mean square value is easily obtained from the exact pairs density [17]: 




(16) 



7T 



Mi : 




I 
1 



Xi(n) (n - r 2 f Xi(r 2 ) dri dr 2 = {2i + 1) 



Xi,i+i{ri) Xi+i,i(r 2 )(ri - r 2 ) 2 dr x dr 2 
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The difference between the interelectronic energy predicted by the Hartree-Fock functional 
(13) and that predicted by the Miiller functional on the true harmonium ground states is 
then given by: 



U 2 — UJ 2 r 

E™(t) - E™(t) = E I 7 * 71 ' - v 7 ^] / Xij(ri) Xjifo) (n - r 2 ) 2 dn dr 2 

i,j=0 J 

oo 

= ~^2[( n2 i -n i )L i + 2{n i n i+1 - v /n^T)M i ] 

i=0 

2ut \l-t 2 1 + t 2 



or equivalently, 



(1 + t) 2 



1 + t 2 (1 + t) 2 

00-/1 



UJ 2 -p 2 



oj + p v Au 

At t — there is no difference between these two values of the energy. It is worth noting 
that there is another point of coincidence, namely t ~ 0.54 or S/k ~ 0.99. Below this value 
the difference is positive, and above it is negative. At t — 1, we find 

Since pf F — p 1 = 2 J2i( n i ~ n i)Xi 7^ for t > 0, this functional does not satisfy the sum 
rule, except when the Hartree-Fock functional is evaluated on a Hartree-Fock state. 



4.3 The Goedecker Umrigar functional 

The Goedecker-Umrigar functional [5] introduces a small variation of Midler's, attempting 
to exclude "orbital self-interaction" . For our closed-shell situation, it is given by: 

p^(r u r 2 ) - pf( ri ,r 2 ) = - n 2 ) Xi (n) Xi (r 2 ). 

i 

This relation implies that for t > this functional violates the sum rule: p\ 7^ 2 j p^ u dr 2 . 
The interelectronic part of the energy difference calculation is given by 

- n 2 )U = [(1 - t 2 r - (1 - t 2 ) 2 t^ U . 

r i 

Hence, the mean value of the interelectronic repulsion predicted by this functional is 



2ut 



1+t 4 
1 -t 4 



1 + t 
1-t 2 



2\2 



(l + t) 2 

The interelectronic energy calculated by means of the Goedecker-Umrigar functional is 
higher than the exact value. At t — 1, the difference diverges. This is unsurprising, given 
that when the coupling is large enough the self-interacting part is almost half of the total 
interelectronic energy; for instance, £^ u (0.8)/£ ee (0.8) ~ 0.44. 
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2-RDM 


A nti svm m etrv 


TTermiticitv 


Sum Rule 


Exact 


YES 


YES 


YES 


Muller 


NO 


YES 


YES 


Hartree-Fock 


YES 


YES 


NO 


GU 


NO 


YES 


NO 


BBC 


NO 


YES 


YES 


CHF 


NO 


YES 


YES 


CGA 


NO 


YES 


YES 



Table 1: Properties fulfilled by the exact two-body functional for two-electron atoms and 
several two-body reduced density approximations [30]. 

4.4 Buijse-Baerends corrected functionals 

A few years after the original Buijse and Baerends' paper [4], some corrections were intro- 
duced, to distinguish between strongly occupied natural orbitals (whose occupation numbers 
are close to 1) and weakly occupied ones (occupation numbers near 0) [31]. The harmonium 
ground state possesses only one strongly occupied orbital, namely xo, whose occupation num- 
ber is n = 1 — t 2 . However, this distinction is lost at high values of the coupling parameter. 
The first corrected functional (BBC1) is given by p 2 BC1 = p¥ + P^i where 



P^ 



ri,r 2 ; 




rimj — 



The second correction (BBC2) modifies BBC1 by adding further terms of the form (^/i 
niiij) Xij{ r i) Xji{ r 2) for distinct strongly coupled orbitals. For the harmonium ground state, 
we may ignore it here; thus we write p 2 BC for p BBC1 . Since both corrections involve only 
off-diagonal terms (r ^ s), these corrected functionals still fulfil the sum rule. 

The functional difference now reads p 2 BC — P2 = p 2 1 and the interelectronic energy 
difference yields 



TTlBBC 771 

-free - free 



M P 2 - 



/ P2 1 r L dri dr 2 = A^2 M 

J i>0 



J2V n i n i+l + = Jl 



\ut A 1 + t 2 



i>0 



+ t) 2 1 - t 2 



1-t 2 



+ 1 



As in the Goedecker-Umrigar functional case, ait — 1 the difference has a divergence. Over 
almost the whole range of t, there is a large error in the energy (see Figure 3). Thus, applied 
to harmonium, these functionals do not reproduce the success found for the homogeneous 
electron gas 
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0,2 0,4 0,6 0,8 1 

t 



Figure 3: The error of the energy value calculated for different functionals. The error 
is defined as the dimensionless [7^ nctlonal (t) _ E cc (t)]/oj, evaluated on the exact one-body 
density matrix for the harmonium ground state. The Miiller functional does not appear here 
since its energy value is exact. 

4.5 CHF and CGA functionals 

Corrected Hartree-Fock (CHF) and Csanyi-Goedecker-Arias (CGA) functionals introduced 
in [6] are improvements of the Hartree-Fock functional. They were designed as tensor prod- 
ucts to get better predictions for the correlation energy in homogeneous electron gases at 
high densities. For a closed shell system, they read 

d,2 HF (u 1 ,u 2 ) = - dx{ux) di(u 2 ) - {rijUj + \Jrii{l - n^n^l - %))xy(«i) Xji( u 2), 

d CGA (u 1 ,u 2 ) = ^di(ux)di(u 2 ) - )^ J (nin j + ^ n<(2 - n i )n j (2 - n^XijM Xa^h)- 

First, note that both functionals satisfy the sum rule: d cnF = d CGA — d\. As regards the 
interelectronic energy, we find that 

£ e c e HF - Ef c = 2 V^^i 1 - V^Ti ~ V(±-^)(l-n i+1 ))M l , 

i 

E GGA -E%= V™^( 2 - V^Ti - V / (2-n i )(2-n i+1 ) ) M t . 

i 

As can be seen in Figures 3 and 4, both functionals show a remarkably good description 
of the energy. At t = and t = 1 the energy is exact. For the CHF functional, the worst 
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GU — ■ — HF 


0,1 
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0,05 
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-0,1 - 




0,1 0,2 0,3 




t 



Figure 4: A zoom of the energy value errors from Figure 3. 



performance occurs around t ~ 0.4 or 5/k ~ 0.96, whose error is (-E£ HF — E^)/ui ~ 0.11; the 
CGA functional is worst at t ~ 0.43 or 5/k ~ 0.97, with an error of (£ e c e HF - Ef e )/u ~ 0.03. 
The estimates of the energy are lower than the correct one; thus they are both overbinding 
for harmonium. 

5 Conclusion 

We have used harmonium as a laboratory to study the performance of some of the one- 
body density functionals proposed to compute the interelectronic repulsion energy in the 
framework of DMFT. We have confirmed the exact value of the energy given by the Miiller 
functional when evaluated on the exact ground state. The functionals which exclude self- 
interacting terms or distinguish between strongly and weakly occupied orbitals display good 
approximation for the energy at small values of the coupling parameter but very poor values 
beyond t ~ 0.3. 

The CHF approximation yields a good description of the interelectronic repulsion, even 
at high values of the parameter t. The performance of the CGA approximation is remarkably 
good, taking into account that it was built explicitly for the electron g clS CclSC The reader 
should keep in mind the violation of some physical constraint or other by each one of the 
examined functionals [33]. 
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